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I. INTRODUCTION 

Studying the production of W pairs at collider exper- 
iments offers a great possibility for tests of the gauge 
sector of the Standard Model, that has been extensively 
investigated by the LEP2 collaborations 0, H 0, Ill- 
Tests in this channel are quite sensitive, because there is a 
destructive interference of two contributions: a f-channel 
contribution, where both W bosons couple to incoming 
fermions, and an s-channel contribution, where the W 
bosons emerge through a triple gauge coupling, cither 
-fW + W~ or ZW + W~ . New physics beyond the Stan- 
dard Model could easily manifest itself, either through 
new particles propagating in the s-channel, like.for in- 
stance, a Z' particle in L-R symmetric models ja, 0,01 m , 
or through anomalous triple gauge couplings, which could 
be loop-induced, mediated b y heavy virtual particles 
running in the loop. In [lol llll Il2j | the most general 
form of an effective Lagrangian for such interactions has 
been developed and discussed. Such tests of anomalous 
trip le g auge couplings have been performed both at LEP2 
El El 01 and at Tevatron, Run I El El El HI 
and at Run II |2l|. Both scenarios could clearly modify 
the total cross section or, at least, lead to different distri- 
butions of the final state particles. In addition, W pairs, 
possibly in association with jets, represent a background 
to a number of relevant other processes, such as the pro- 
duction of top quarks, the production of a Higgs boson 
with a mass above roughly f 35 GeV, or the production of 
supersymmetric particles, such as charginos or neutrali- 

nos nim. 

Accordingly, there are a number of calculations and pro- 
grams dealing with this process. At next-to-leading order 
(NLO) in the strong couplin g co nst ant. W pair produc- 
tion has been calculated by |24 EM |26| . In addition, a 
number of programs have been made available, allowing 
the user to implement phase space cuts and to generate 
single events. First of all, there are fixed order calcu- 
lations. At leading order (LO), i.e. at tree- level, they 
are usually performed through automated tools, called 
matrix element or parton level generators. Examples 
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for such programs include COMPH EP [27j . GRACE/GROPPA 
[HHIl, MADGRAPH/MADEVENT [30l l3lf 7~ALPGEN 32], and 
AMEGIC++ 33). At NLO, the program MCFM [34J provides 
cross sections and distributions for this process. 
Apart from such fixed order calculations, multipurpose 
eve nt g enerators such as PYTHIA jH, or HERWIG 
|3?1 l38j | play a major role in the experimental analyses 
of collider experiments. They proved to be extremely 
successful in describing global features of such processes, 
like, for instance, the transverse momenta or rapidity dis- 
tributions of the bosons. They are usually based on exact 
tree-level matrix elements for the production and decay 
of the boson pair, supplemented with a parton shower. 
The latter takes proper care of multiple parton emission 
and resums the corresponding leading and some of the 
subleading Sudakov logarithms. 

In view of the need for increasing precision, recently two 
approaches have been developed that incorporate higher 
order corrections into the framework of multipurpose 
event generators. The first one, called MCONLO, provides a 
method to consistently match NLO calculations for spe- 
cific processes with the parton shower |39|, |40j . The idea 
of this approach is to organize the counter-terms neces- 
sary to cancel real and virtual infrared divergencies in 
such a way that the first emission of the parton shower 
is recovered. Of course, this method depends to some 
extent on the details of the parton shower, and it has 
some residual dependence on the process in question. So 
far, MCONLO has been implemented in conjunction with 
HERWIG |4f| for the following processes: production of W 
and Z bosons, or pairs of these bosons (3!|, production 
of the Higgs boson, production of heavy quarks |4fJ . 
An alternative approach is to consistently combine tree- 
level matrix elements for different multiplicities of addi- 
tional jets and to merge them with the parton shower. 
This approach has been presented for the first time for 
the case of e + e~ annihilations into iets [421: later it 
has been extended to hadronic collisions [43| and it has 
been reformulated to a merging procedure with a dipole 
shower in |44| . The idea underlying this method is to 
separate the kinematical range of parton emission by a 
fc_i_-algorithm 0, EE into a regime of jet produc- 
tion, covered by the appropriate matrix elements, and a 
regime of jet evolution, covered by the respective shower. 
Then, the matrix elements are reweighted through Su- 
dakov form factors and hard emissions in the parton 
shower leading to a jet are vetoed such that there is only 
a residual dependence on the jet resolution cut. This 
method is one of the cornerstones of the new event gen- 
erator SHERPA 48] ; it has been validated for the cases of 
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e+e annihilations into jets 0,H3 and for the produc- 
tion of single vector bosons at the Fermilab Tevatron 
and the CERN LHC 

In this publication this series of studies will be contin- 
ued with an investigation of W pair production at the 
Fermilab Tevatron, Run II, where both W bosons decay 
lcptonically, i.e. pp — > W + W~ + X — ► e + n~v e v^ + X . 
Input parameters used throughout this publication and 
the specifics, how the SHERPA runs have been obtained, 
are listed in the appendix, see Apps. El and 10 After 
some consistency - including scale variation - checks of 
the merging algorithm in Sec. [HJ results obtained with 
SHERPA will be confronted with those from an NLO cal- 
culation provided by MCFM, cf. Sec. IIIII Then, in Sec. 
II VI some exemplary results of SHERPA are compared with 
those obtained from other event generators, in particular 
with those from PYTHIA and MCONLO. A summary closes 
this publication. 



II. CONSISTENCY CHECKS 

In this section some sanity checks of the merging algo- 
rithm for the case of W pair production are presented. 
For this, first, the dependence of different observables on 
the key parameters of the merging procedure, namely the 
internal matrix-element parton-shower separation scale 
Q cu t and the highest multiplicity n max of included tree- 
level matrix elements, is examined. Secondly, the sensi- 
tivity of the results with respect to changes in the renor- 
malization scale /j,r and the factorization scale [If will 
be discussed. 

All distributions shown in this section are inclusive re- 
sults at the hadron level, where restrictive jet and lepton 
cuts have been applied, for details on the cuts cf. App. 
lO In all cases, the distributions are normalized to one 
using the respective total cross section as delivered by 
the merging algorithm. 



Impact of the phase space separation cut 

First of all, the impact of varying the jet resolution cut 
Q cu t is studied. SHERPA results have been obtained with 
an inclusive 2jet production sample, i.e. tree-level ma- 
trix elements up to two additional QCD emissions have 
been combined and merged with the parton shower. In 
all figures presented here the black solid line shows the 
total inclusive result as obtained by SHERPA for the re- 
spective resolution cut Qcut- The reference curve drawn 
as a black dashed line has been obtained as the mean 
of five different runs, where the resolution cut has been 
gradually increased, Q cu t = 10, 15, 30, 50 and 80 GeV. 
The coloured curves represent the contributions stem- 
ming from the different matrix- clement final-state multi- 
plicities. Results are shown for three different resolution 
cuts, namely Q cut = 15, 30 and 80 GeV. It should be 
noted that the change of the rate predicted by the merg- 
ing procedure under Q cu t variation has been found to be 



very small, although it is a leading order prediction only. 
Nevertheless, by varying the separation cut between 10 
and 80 GeV, the deviation of the total rate amounts to 
2.4% only. 

As a first result, consider the pt distribution of the W + 
boson, presented in Fig. ^ The distributions become 
slightly softer for increasing cuts. However, this observ- 
able is very stable under variation of Q C ut with maximal 
deviations on the ±5% level only. The shape of the W + 
boson's pt is already described at LO (using a parton 
shower only). As it can be seen from the figure, this LO 
dominance is nicely kept by the SHERPA approach under 
Q cu t variation. There the ljet (green line) and 2jet (blue 
line) contributions are reasonably - for the 80 GeV run, 
even strongly - suppressed with respect to the leading 
contribution. 

In Fig. |21 the transverse momentum spectrum of the 
W + W~ system is depicted. Here, deviations show up, 
but they do not amount to more than ±20%. Thus, the 
QCD radiation pattern depends only mildly on Q cu t (in- 
dicated by a vertical dashed-dotted line), which at the 
same time has been varied by nearly one order of mag- 
nitude. For Q cut = 15 GeV the matrix element domain 
is enhanced with respect to the reference resulting in a 
harder px tail. In contrast by using Q cut = 80 GeV the 
hard tail of the diboson transverse momentum is under- 
estimated with respect to the reference, since the parton 
shower attached only to the lowest order matrix element 
starts to fail in the description of high-p^ QCD radiation 
at pt ~ 30 GeV. At Q cut = 80 GeV a smooth transition 
is required. The higher order matrix elements then stop 
the decrease in the px prediction. 

In previous publications it turned out that differential 
jet rates most accurately probe the merging algorithm, 
since they most suitably reflect the interplay of the ma- 
trix elements and the parton shower in describing QCD 
radiation below and above the jet resolution cut. Results 
obtained with the Run II fc^-algorithm using R = 1 are 
shown for the 1 — > 0, 2 — ► 1 and 3^2 transition in the 
left, middle and right panels of Fig. |3 respectively. The 
value for the internal cut increases from Qcut = 15 GeV 
(top) to Qcut = 80 GeV (bottom). Compared with the 
Pjf w spectra, similar characteristics of deviations from 
the reference curve appear. However, here, they are mod- 
erately larger reaching up to ±30%. The dashed dotted 
vertical line again marks the position of Qcut, which also 
pictures the separation of the njet from the n + ljet con- 
tribution. Small holes visible around the respective sep- 
aration cuts are due to a mismatch of matrix element 
and parton shower kinematics. For Q cut = 80 GeV these 
holes are much more pronounced, reflecting the failure of 
the parton shower in filling the hard px emission phase 
space appropriately. 

Taken together, the deviations found are very moderate; 
however, in certain phase space regions they may reach 
up to 30%. This is satisfactory, since the merging algo- 
rithm guarantees Qcut independence on the leading log- 
arithmic accuracy only. The residual dependence of the 
results on Qcut may be exploited to tune the perturbative 
part of the Monte Carlo event generator. 



Singly resonant diagrams contributing to the parton level pro- 
cesses of pp — > e+ \i~ VcDp, + X have been included. 
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FIG. 1: The pr distribution of the W boson and its depen- 
dence on Qcut, chosen to be 15, 30 and 80 GeV (from top to 
bottom). The black solid line shows the SHERPA prediction 
obtained with n max = 2, the black dashed one is the reference 
obtained as the mean of different Q C ut runs and the coloured 
lines indicate the different multiplicity contributions. The 
lower part of the plots exhibits the normalized difference of 
the prediction with respect to the reference. Cuts and input 
parameters are specified in the appendices. 



FIG. 2: The pr distribution of the W W~ system under 
merging scale variation. The cut indicated through a vertical 
dashed-dotted line has been chosen as Q cu t = 15, 30 and 80 
GeV (from top to bottom). The black solid line shows the 
SHERPA prediction obtained with n max = 2, the black dashed 
one is the reference obtained as the mean of different Q C ut 
runs and the coloured lines indicate the different multiplicity 
contributions. The lower part of the plots exhibits the normal- 
ized difference of the prediction with respect to the reference. 
Cuts and input parameters are specified in the appendices. 



4 




FIG. 3: Differential 1 — > jet rate Qi, 2 — > 1 jet rate Q2 and 3 — > 2 jet rate Q3 (left to right) for the SHERPA n max = 2 
configuration. The cut has been chosen to be 15, 30 and 80 GeV (from top to bottom). The black solid line shows the total 
result, the black dashed one is the reference obtained as the mean of different Q C ut runs and the coloured lines indicate the 
different multiplicity contributions. The vertical dashed dotted line indicates the separation cut position. The lower part in 
all plots pictures the normalized difference of the corresponding prediction with respect to the reference. For input parameters 
and cuts, see Adds . lAl and IU1 



Impact of the maximal number of included matrix 
elements 

The approach of varying the maximal jet number n max 
can be exploited to further scrutinize the merging pro- 
cedure. In all cases considered here, Qcut has been fixed 
to Qcut = 15 GeV. This maximizes the impact of higher 
order matrix elements. In spite of this, for very inclusive 
observables, the rates differ very mildly, the change is less 
than 2%. 

In Fig. 01 once more the transverse momentum distribu- 
tion of the W + gauge boson is presented, illustrating that 
the treatment of the hig hest multiplicity matrix elements 
(for more details cf. [5(J, |51| ) completely compensates for 



the missing 2jet matrix element in the n max = 1 case. 
The behaviour is almost unaltered when changing from 
the n max = 1 to the n max = 2 prediction (cf. the right 
panel). In contrast, n max = yields a considerably softer 
distribution (cf. the left panel). 

Lepton pt spectra show similar characteristics like the 
W + distribution. However, there are a number of ob- 
servables, which turned out to be rather stable under the 
variation of n max , such as the pseudo-rapidity spectra of 
the W + boson, the positron and muon or correlations 
between the leptons, e.g. the A<fi or Ai? distribution. 
In these cases, deviations turn out to be smaller than 
±5% in total, i.e. when considering the change between 
the pure shower and the inclusive 3jet production per- 
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FIG. 4: The pr distribution of the W + boson in dependence on the variation of the maximal jet number. The comparison is 
to a (black dashed) reference curve obtained with n£; ax = n max — 1. The cut has been chosen to be 15 GeV. In both plots the 
black solid line shows the total result obtained with SHERPA. The coloured lines indicate the different multiplicity contributions. 
The lower part in both plots visualizes the normalized difference of the corresponding prediction with respect to the reference. 
For input parameters and cuts, see Adds . lAl and IU1 
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FIG. 5: The Ht distribution and its dependence on the varia- 
tion of the maximal jet number. The separation cut has been 
chosen to be 15 GeV. The green solid line shows the SHERPA 
prediction obtained with n max = 1, the lighter dashed and the 
grey dotted one stand for the n max = 2 and n max = 3 predic- 
tion, respectively; the darkgreen dashed-dotted curve pictures 
the pure shower performance of SHERPA starting off with the 
lowest order matrix element. The lower part of the plot shows 
the normalized differences with respect to the n max = 1 case. 
For input parameters and cuts, see Adds . lAl and IU1 



formance of SHERPA. Even the pseudo-rapidity spectra of 
the resolved jets are rather unaffected. 
In contrast, three more observables are presented show- 
ing a sizeable (< ±30%) or even strong (ps ±100%) de- 
pendence on the variation of the maximal jet number, 
namely the Ht distribution depicted in Fig. and the 
inclusive pt spectra of associated jets exhibited in Fig. 
The upper and lower panel of Fig. shows the spectra 
of the hardest and the second hardest jet, respectively. 



Owing to the nature of these three observables to be sen- 
sitive on extra jet emissions, predictions - as expected 
- become harder with the increase of n max . However, a 
stabilization of the predictions is clearly found with the 
inclusion of more higher order matrix elements describing 
real QCD emissions. 



Effects of renormalization and factorization scale 
variations 

In the following the impact of renormalization and fac- 
torization scale variations is discussed. For the SHERPA 
merging approach, this variation (also cf. H2) * s P er ~ 
formed by multiplying all scales with a constant factor 
in all coupling constants and PDFs, which are relevant 
for the matrix element evaluation, the Sudakov weights 
and for the parton shower evolution. 
For this study, the SHERPA samples are produced with 
n max = 1 and Qcut = 15 GeV. In all figures the 
green solid line represents SHERPA's default scale choices, 
whereas the black dashed and the black dotted curve 
show the outcome for scale multiplications by 0.5 and 2.0, 
respectively. The total rate as provided by the merging 
algorithm is again remarkably stable, varying with re- 
spect to the default only by ±4.2%, thereby increasing 
for smaller scales. 

The transverse momentum distribution of the W + boson 
is investigated in Fig. Scale variations slightly distort 
the shape, shifting it towards harder px for smaller scales 
and vice versa. The effect is more pronounced in the Ht 
distribution, shown in Fig. |H1 and in the transverse mo- 
mentum distribution of the diboson system, depicted in 
Fig. [5] However, the deviations maximally found reach 
up to ±30%. In contrast to the findings stated so far, 
jet transverse momentum spectra do not feature shape 
distortions under scale variations. 

The pattern found from these investigations can be ex- 
plained as follows. The single matrix element contribu- 
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FIG. 6: SHERPA predictions of the inclusive pr of the asso- 
ciated jets considered in dependence on the variation of the 
maximal jet number. The spectra of the hardest and the sec- 
ond hardest jet are depicted in the upper and the lower panel, 
respectively. The jet resolution cut has been taken to be 15 
GeV. The green solid line shows the result of the n max = 1 
sample, the brighter dashed and the grey dotted one stand 
for the n max = 2 and n max = 3 sample, respectively; the 
darkgreen dashed-dotted curve depicts the pure shower per- 
formance. The lower part of the plot shows the normalized 
differences with respect to the n max = 1 case. For the jet def- 
inition, the Run II fc^-algorithm with R = 0.7 and > 15 
GeV has been used. For more details, see Adds. 1X1 and ICl 

tions - here the Ojet and ljet contribution - have their 
own rate and shape dependencies under scale variations. 
In their interplay these differences transfer to changing 
the admixture of the single contributions. Hence, shape 
modifications can appear as soon as different phase space 
regions are dominated by a single contribution. This also 
explains the behaviour found for jet prs. In the case 
studied here, they are solely described by the ljet ma- 
trix element with the parton shower attached, thus, their 
different rates cancel out due to normalization and their 
shapes are not affected. 

Taken together, the dependencies found here, together 
with the ones on Qcut and n max , yield an estimate for 
the uncertainty related to the SHERPA predictions. 
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FIG. 7: The pr distribution of the W + boson under scale 
variations. All predictions stem from SHERPA with n max = 1 
and Qcut = 15 GeV. The green solid line shows the prediction 
under default scale choices for the merging procedure. For 
the black dashed and the black dotted curve, all scales for 
the coupling constants and PDFs have been multiplied by 
0.5 and 2.0, respectively. The lower part of the plot presents 
the normalized differences with respect to the default choice. 
Input parameters and cuts are given in Adds . lAl and IgI 
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FIG. 8: The Ht distribution and its dependence on the vari- 
ation of fiR and fiF in the merging prescription. Fixing 
i max = 1 and Qcut = 15 GeV, the green solid line shows the 
prediction under default scale choices. For the black dashed 
and the black dotted curve, all scales for the coupling con- 
stants and PDFs have been multiplied by 0.5 and 2.0, respec- 
tively. The lower part of the plot presents the normalized 
differences with respect to the default choice. Input parame- 
ters and analysis cuts are given in Adds. lAl and IU1 



III. SHERPA COMPARISON WITH MCFM 

In this section, the focus shifts from internal sanity checks 
to comparisons with a full NLO calculation. For this, 
the MCFM program [34| has been used. In both, MCFM and 
SHERPA the CKM matrix has been taken diagonal, and 
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FIG. 9: The pr distribution of the W pair under variation of 
Hr and hf- Fixing n max = 1 and Q cu t = 15 GeV, the green 
solid line shows the prediction under default scale choices. 
The black dashed and the black dotted curve is generated 
when all scales used for the coupling constants and PDFs have 
been multiplied by 0.5 and 2.0, respectively. The lower part 
of the plot presents the normalized differences with respect to 
the default choice. Input parameters (including a primordial 
k± smearing) and cuts are given in Adds. 1X1 and Icl 



no b quarks are considered in the partonic initial state 
of the hard process. If not stated otherwise, in MCFM the 
renormalization and factorization scale have been chosen 
as fiR = i-iF = Mw, according to the choice made in [34| . 
For more details on the input parameters and setups, 
see Apps. 1X1 andlEl In the following the results of MCFM 
are confronted with those of SHERPA (using Q cu t = 15 
GeV) obtained at the parton shower level. Furthermore, 
for this analysis, realistic experimental cuts (cf. App. IC)l 
have been applied and all distributions have been nor- 
malized to one. 

First the Ht distribution, depicted in Fig. is con- 
sidered. Clearly, higher order corrections affect the Ht 
shape. This is due to two reasons. First of all, the addi- 
tional QCD radiation may manifest itself as jet(s), which 
thus contribute to Ht- Otherwise the additional partons 
still form a system against which the W pair may recoil. 
Quantitatively, the inclusion of NLO results in a shift 
of the Ht distribution at harder values by up to 20%; 
in SHERPA this trend is amplified by roughly the same 
amount. The differences between MCFM and SHERPA, how- 
ever, are due to the different scale choices in both codes. 
In MCFM all scales have been fixed to /i = Mw, whereas, 
forced by the merging procedure, in SHERPA the scales 
are set dynamically. In view of the scale variation results 
discussed in the previous section for Ht (cf. Fig. [8} devi- 
ations of this magnitude owing to different scale choices 
are possible. 

The impact of scale variations on the shape of the same 
observable is quantified in Fig. ^2 This time, however, 
the SHERPA result with n max = 1 is compared to NLO re- 
sults obtained from MCFM with scale choices in the range 
fJ-B. = fJ-F = Mw ■ ■ -4:M\y and with a LO result taken 
at = fip = 2Mw- Obviously, the smaller choice of 
scale results in the MCFM outcome to be closer to the one 
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FIG. 10: Normalized H T distribution. SHERPA results are 
shown for n max = 1 (green solid line) and n max = 2 (green 
dashed line) and compared to the QCD NLO result of MCFM 
(black solid line). The LO result with the same scale choice is 
depicted as a thin black dashed line. A difference plot with the 
MCFM NLO prediction as reference is given within the figure. 



of SHERPA. As expected, in comparison to the scale vari- 
ation results found for SHERPA, the shape uncertainties 
of the full NLO prediction due to varying the scales are 
smaller. 

In Fig. 1121 Ht is depicted again, this time for the case 
of exclusive pp — > e + [i^v^D^ production. There, the real 
part of the NLO correction in MCFM is constrained such 
that it does not produce an extra jet (for jet definition, 
see App. [UJ|. In SHERPA the Ojet matrix element with 
the parton shower attached is considered exclusively, i.e. 
the parton shower is now forced not to produce any jet 
at all. In this case, the higher order corrections lead to 
a softer Ht distribution compared to the leading order 
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FIG. 11: Normalized Ht distribution. Here both, the renor- 
malization and factorization scale of the NLO calculation have 
been varied in the range {ir = [If = Mw ■ • ■ 4 Mw , indicated 
by the shaded area. These MCFM results are compared with the 
leading order result at [lr = /if = 2 Mw (thin black dashed 
line) and with the result of SHERPA where n max = 1 (green 
solid line). The lower part of the plot shows the normalized 
differences with respect to the SHERPA result. 
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FIG. 12: Normalized Ht distribution in exclusive pp — > 
e + pL~v e Vfj, production. The SHERPA result (green solid line) 
is obtained with n max = and a parton shower constrained 
not to produce any extra jets. This result is compared with 
the MCFM result at NLO in a a (black solid line) and with the 
LO result (thin black dashed line). The latter two are taken 
for the default scale choices. The lower part of the plot shows 
the normalized differences with respect to the NLO result ob- 
tained from MCFM. 



FIG. 14: Normalized Af&^^y distribution of the W boson sys- 
tem. The MCFM result (black line) is contrasted with results 
from SHERPA at the parton shower level with n max = (dark- 
green dotted line), n max = 1 (green solid line) and n ma x = 2 
(green dashed line) . Again a primordial k± smearing has been 
used. Additionally, the blue dashed curve represents a predic- 
tion obtained with MC@NL0. The lower part of the plot shows 
the normalized differences with respect to the result of MCFM. 
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FIG. 13: Normalized pr distribution of the W pair. The 
MCFM result at fin = /if = Mw (black line) is contrasted with 
the predictions made by SHERPA both at the matrix element 
level (blue line) and at the parton shower level with n max = 
(darkgreen dotted line) and n max = 1 (green solid line). A 
primordial k± smearing has been used to obtain the SHERPA 
shower results. 
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FIG. 15: Normalized transverse momentum distribution of 
the W + boson. The results of SHERPA for n max = 1 (green 
solid line) and for n max = 2 (green dashed line) are compared 
with the QCD NLO result obtained by MCFM (black solid line) 
and with the LO result (thin black dashed line) for the default 
scale choices, i.e. {ir = /if = Mw- Within the plot the 
normalized differences with respect to the NLO result of MCFM 
are given. 



prediction, and the results of MCFM and SHERPA show the 
same deviations as before (cf. Fig.EJ- 
The effect of QCD radiation is best observed in the pr 
distribution of the W pair, depicted in Fig. ED Clearly, 
without any radiation, the px of the W pair is exactly 
zero, and only the emission of partons leads to a recoil 
of the diboson system. In the NLO calculation of MCFM, 
however, the spectrum is therefore described at lowest 
order, in this particular case taken at /ir = fiF = Myy- 
In contrast, in the SHERPA matrix element result, sub- 
jected to the explicit jet cut, Sudakov form factors and ot s 
reweighting are applied with a variable scale choice, ex- 
plaining the differences between the two matrix-element 



type results in this figure. Contrasting this with the par- 
ton shower approach, it is clear that parton emission 
through the shower alone is not sufficient to generate 
sizeable px of the W pair in the hard region. For this, 
the corresponding matrix element has to be employed, 
leading to a very good agreement with the MCFM outcome 
in the high-pT tail of the distribution. In the soft regime 
the result of the bare MCFM matrix element is unphysical. 
Due to the cascade emission of soft and collinear partons, 
SHERPA accounts for resummation effects, which clearly 
yield the depopulation of the softest-pr region. 
Another way to look at the effects of QCD radiation is to 
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FIG. 16: Normalized transverse momentum distribution of 
the e + produced in the decay of the W + . The results of 
SHERPA for ^max - — 1 (green solid line) and for 
(green dashed line) are confronted with the QCD NLO result 
obtained by MCFM (black solid line) and with the LO result 
(thin black dashed line). For the latter two, the scales are 
again fixed according to the default choices, i.e. /j,r = [if = 
Mw- Within the plot the normalized differences with respect 
to the NLO result of MCFM are shown. 



FIG. 18: Normalized AR distribution between the two 
charged leptons, the positron and the muon, emerging from 
the W decays. SHERPA results for n max = 1 (green solid line) 
and n max = 2 (green dashed line) are compared to those pre- 
dicted by MCFM (black solid line). The LO result with the 
same scale choice, is shown as a black dashed line. The lower 
part of the plot shows the normalized differences with respect 
to the NLO result of MCFM. 
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FIG. 17: Normalized r\ distribution of the W boson. The 
SHERPA results for n max = 1 (green solid line) and n max = 2 
(green dashed line) are confronted with those of MCFM (black 
solid line) and with the LO result (thin black dashed line). 
Again, in the latter two the scales are chosen as /j,r = (j,f = 
Mw- The normalized differences with respect to the NLO 
result of MCFM are also shown. 



consider the relative angle between the two W bosons 2 , 
see Fig. El Of course, when they decay into leptons plus 
neutrinos this is not an experimental observable, on the 
generator level, however, it is very nice to visualize the 
effect of QCD radiation in this way. Without any QCD 
radiation, the two Ws would be oriented back-to-back, at 
A$^ w — ir. Including QCD radiation, this washes out, 
as depicted in the figure. Again, resummation effects al- 



ter the result of the matrix element alone by decreasing 
the amount of softest radiation, this time corresponding 
to the back-to-back region around A<I>* ss ir. The effect 
of high-px radiation can be clearly seen for small A$* by 
comparing the different n max predictions of SHERPA. The 
larger n max is chosen, the harder the prediction for small 
A<T>*. On the other hand to better value the influence 
of the parton shower a prediction made by MC@NL0 (see 
Apd. lB|l has been included. For a wide region of A<£>*, it 
well agrees with the SHERPA result for n max = 1. 
Figs. Upl and 1161 exhibit the transverse momentum distri- 
butions of the W + and of the e + produced in its decay, 
respectively. Only mild deviations less than 10% between 
MCFM and SHERPA are found, which again can be traced 
back to different scale choices in both approaches. These 
differences recur as and, therefore, explain part of the 
deviations found in the Ht spectrum, cf. Fig. ^] As ex- 
pected, the inclusion of the 2jet contribution in SHERPA 
gives no further alterations of the n max = 1 result. Of 
course, the different radiation patterns also have some 
minor effects on the 77 distribution of the W + depicted 
in Fig. El In the AR efl distribution presented in Fig. 
IThI the NLO result of MCFM and the parton shower level 
results of SHERPA are in nearly perfect agreement with 
each other. Higher order effects tend to change the shape 
of the LO prediction with respect to the NLO one by 
roughly 10%. The interesting observation here is that 
this change is seemingly not related to the transverse 
hardness of a jet system against which the W pair re- 
coils. This gives rise to the assumption that the change 
with respect to the LO result is due to some altered spin 
structure in the 2 — > 5 matrix clement. 



IV. COMPARISON WITH OTHER EVENT 
GENERATORS 



2 The angle is measured in the frame, where the W + W~ system 
rests at the beam axis, i.e. the diboson system is corrected on its 
initial z boost. 



In this section a comparison of SHERPA with other hadron 
level event generators, in particular PYTHIA and MC0NL0 
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will be discussed. Details on how their respective samples 
have been produced can be found in the Adds. 1X1 and FBI 
The SHERPA samples have been generated with n max = 1 
and Q cu t = 15 GeV. The comparison is again on inclusive 
distributions - normalized to one - under the influence 
of realistic experimental cuts, for details see App. IU1 



Comparison of the QCD activity 

As before, the starting point is the discussion of the ra- 
diation activity predicted by the various codes. In Fig. 
1191 results for the Ht observable obtained from PYTHIA, 
MC@NL0 and SHERPA are displayed. The predictions of the 
former two codes nicely agree with each other. Similar 
to the SHERPA MCFM comparison, SHERPA again predicts a 
slightly harder spectrum, with relative deviations of up 
to 20%. 

Closer inspection of the reason for the differences in the 
Ht spectrum reveals that the agreement of PYTHIA and 
MCQNLO is presumably a little bit accidental. A first hint 
into that direction can be read off Fig. 1201 where the 
Pt spectrum of the W pair is displayed. In the region 
of low p T (up to 100 GeV), the results of MCONLO and 
SHERPA are in fairly good agreement 3 , and sizeable dif- 
ferences larger than 10% appear only for pr > 100 GeV. 
In contrast, the PYTHIA result for this observable shows 
a significant enhancement of the low-px region and stays 
well below the other predictions for pt > 10 GeV. This 
comparison of the three differential cross sections clearly 
underlines that the three codes differ in their description 
of the QCD emissions. 
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FIG. 19: Normalized Ht distribution obtained from PYTHIA 
(red dotted line) , MCONLO (blue dashed line) and SHERPA (green 
solid line). For the generation of the SHERPA sample, n max = 1 
and Qcut = 15 GeV have been chosen. The lower part of the 
plot exhibits the normalized differences with respect to the 
SHERPA prediction. Input parameters and the employed cuts 
are specified in the Adds. 1X1 and IU1 



3 Apart from the very soft region, where the difference is due to 
parton shower cutoff effects in HERWIG. 



Fig. |^ depicts the norm of the scalar difference of the 
transverse momenta of the W + and W~ gauge boson, 
\pr — Pt I • This observable is sensitive to higher or- 
der effects, since at LO it merely has a delta peak at 
Pt = GeV. Again, the hardest prediction is delivered 
by SHERPA with n max = 1, results from MCONLO, PYTHIA, 
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FIG. 20: Normalized pr distribution of the W + W~ system. 
Results from PYTHIA (red dotted line), MCONLO (blue dashed 
line) and SHERPA (green solid line) are compared. For the 
generation of the latter, n max = 1 and Q cu t = 15 GeV have 
been chosen. The lower part of the plot presents the normal- 
ized differences with respect to the SHERPA prediction. Input 
parameters (including a primordial k± smearing) and the em- 
ployed cuts are specified in the Adds. 1X1 and IU1 




FIG. 21: Difference of the scalar transverse momenta of the 
two W bosons, \pr — Pt I ■ The predictions compared are: 
PYTHIA given as a red dotted curve, MCONLO depicted by the 
blue dashed line and SHERPA in inclusive ljet production at 
Qcut = 15 GeV drawn as a green solid line as well as SHERPA 
in pure shower performance shown as a darkgreen solid line. 
The lower part of the plot shows the normalized differences 
with respect to the SHERPA prediction with n max = 1. In- 
put parameters and the employed cuts are summarized in the 
Adds. 1X1 and IU1 
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FIG. 22: Transverse momentum distributions of the associ- 
ated jets, in the upper panel, the inclusive pr of the hardest 
jet is depicted, whereas in the lower panel that one of the sec- 
ond hardest jet is displayed. Again, results from PYTHIA are 
given by the red dotted lines, MC@NL0 results are represented as 
blue dashed lines and SHERPA results are the green solid lines. 
For the generation of the latter, n max = 1 and Q cu t = 15 GeV 
have been used. The lightgreen solid line in the lower panel 
corresponds to the SHERPA result obtained with n max = 2. 
The lower part of both plots shows the normalized differences 
with respect to the SHERPA n max = 1 performance. The in- 
put parameters and the employed cuts are summarized in the 
Apps.|X]and|0 



and the pure shower performance of SHERPA are increas- 
ingly softer. For |Apr| > 60 GeV, this observable seems 
to depend more and more on the quality of modelling the 
hardest emission, which is intrinsically better described 
by MCONL0 and by SHERPA with n max = 1. The fact that 
the PYTHIA shower performs better than the pure SHERPA 
shower for high px differences can be traced back to the 
choice of starting scale for the shower evolution, which is 
either s pP (PYTHIA) or s W w (SHERPA). 
In fact, differences appear in the px distributions of the 
hardest two jets, see Fig.|221 The upper part of this fig- 
ure depicts the transverse momentum spectrum of the 
hardest jet. Surprisingly, although MC0NL0 contains a 
matrix element for the emission of an extra jet, its px 



distribution is considerably softer (by up to 40%) than 
the result of SHERPA generated with n max = 1. This 
trend is greatly amplified when going to the spectrum of 
the second hardest jet. There, clear shape differences of 
the order of a factor 2 between the SHERPA ljet sample 
and MC0NL0 show up for px « 180 GeV. The surprise 
according to this figure is that PYTHIA and SHERPA using 
n max = 1 almost agree on the px distribution of the sec- 
ond jet, although they were different for the hardest jet. 
At that point it should be noted that the second jet in 
both cases, PYTHIA and SHERPA with n max = 1, is pro- 
duced by the parton shower only. Given the drastically 
larger shower start scale of PYTHIA, it seems plausible 
to achieve to some extent a compensation for the intrin- 
sic parton shower deficiencies in filling the hard emis- 
sion phase space 4 . However, in the very moment, SHERPA 
events are generated with appropriate matrix elements, 
i.e. with n max = 2, this distribution is dramatically dif- 
ferent for the three codes with deviations larger than a 
factor 2 for p T « 120 GeV. 

Taken together, these findings hint that the three codes 
differ in their modelling of the QCD activity, especially 
in those of the hardest QCD emission. For MC@NL0 and 
SHERPA the latter can be traced back to the different 
ansatz in including the matrix element for this emission, 
where again different scale choices may trigger effects on 
the 20% level. 



Comparison of lepton observables 

Finally, the leptons in the final state as described by the 
three event generators PYTHIA, MC0NL0 and SHERPA will 
be investigated. There, some significant differences ap- 
pear between SHERPA and PYTHIA on the one hand, and 
MC0NL0 on the other hand. These differences are due 
to the fact that at the moment spin correlations of the 
W decay products are not implemented in MC0NL0 5 . To 
validate that effects are indeed due to the lack of spin 
correlations, SHERPA samples have been prepared, where 
these correlations are artificially switched off. Further- 
more, in order to quantify these effects without any bias, 
results have been obtained without the application of any 
lepton and jet cuts. 

The impact of the lack of spin correlations already be- 
comes visible in one-particle observables, such as the px 
or the r\ spectrum of the positron produced in the W + 
decay. These are shown in Figs. 1231 and 1241 respectively. 
Confronting the two methods with each other, which cor- 
rectly respect spin correlations, for the transverse mo- 
mentum distribution of the e + , the following pattern is 
revealed. Due to the consistent inclusion of higher order 
tree-level matrix elements, the SHERPA n max = 1 setup 
produces a considerably harder spectrum than PYTHIA. In 
contrast, the distributions with no spin correlations both 
result in an even harder high-p^ tail. They agree quite 
well up to px = 60 GeV, hence, this coincidence may 



4 PYTHIA's ability to account for harder second jets with respect to 
MCONL0 is a hint for the similarity of their Ht predictions. 

J This situation is currently being cured by the authors of MCONL0 
who prepare a new version of their code including spin correla- 
tions |53l . 
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be assigned to the lack of spin correlations in the gauge 
boson decays. Above that region, the MCONLO spectrum 
again becomes softer with respect to the SHERPA predic- 




100 120 140 

p T (e*)[GeV] 



FIG. 23: Normalized pr spectrum of the positron. Results of 
PYTHIA (red dotted line) and SHERPA (green solid line) includ- 
ing spin correlations are confronted with those obtained from 
MCSNLO (blue dashed line) and with results from SHERPA, where 
spin correlations have been switched off (green dashed line). 
All predictions are generated without the use of cuts. The 
vertical dashed-dotted line is added to indicate the position 
of the usually employed lepton pr cut. For input parameters, 
see App.[X] The lower part of the plot shows the normalized 
differences with respect to the SHERPA prediction including 
spin correlations. 
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FIG. 24: Normalized r\ spectrum of the positron. Results 
of PYTHIA (red dotted line) and SHERPA (green solid line) in- 
cluding spin correlations are compared with those obtained 
from MC@NL0 (blue dashed line) and with results from SHERPA, 
where spin correlations have been switched off (green dashed 
line). All predictions are generated without any restriction. 
The vertical dashed-dotted lines are added to indicate the 
position of the usually employed lepton r\ cuts. For input pa- 
rameters, see App. ^ The lower part of the plot shows the 
normalized differences with respect to the SHERPA prediction 
including spin correlations. 



tion where the spin correlations have been eliminated. 
The fact that all four predicted distributions alter in their 
shape is not solely triggered by the different spin correla- 
tion treatments, again, the different descriptions of QCD 
radiation clearly contribute to the deviations found. 
In contrast, a simpler pattern is found for the aforemen- 
tioned rj distribution of the e + . The results of PYTHIA 
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FIG. 25: Normalized A<^ eM distribution. Results of PYTHIA 
(red dotted line) and SHERPA (green solid line) including spin 
correlations are compared with those obtained from MC@NL0 
(blue dashed line) and with results from SHERPA, where spin 
correlations have been switched off (green dashed line). All 
predictions are obtained without the use of cuts. For input 
parameters, see App.^ The lower part of the plot shows the 
normalized differences with respect to the SHERPA prediction 
including spin correlations. 
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FIG. 26: Normalized A7? eM distribution. Results of PYTHIA 
(red dotted line) and SHERPA (green solid line) including spin 
correlations are compared with those obtained from MC0NL0 
(blue dashed line) and with results from SHERPA, where spin 
correlations have been eliminated (green dashed line). All 
predictions are obtained without the use of cuts. For input 
parameters, see App.^ The lower part of the plot shows the 
normalized differences with respect to the SHERPA prediction 
including spin correlations. 
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FIG. 27: In the left and middle panel the pr spectrum of the W before and after the application of cuts is depicted, respectively. 
The right panel exhibits the r\ distribution of the W + under the influence of these cuts. The predictions compared are: PYTHIA 
(red dotted line), SHERPA (green solid line), MC@NL0 (blue dashed line) and SHERPA without correlations in the boson decays 
(green dashed line). For input parameters, see App. [X] The lower part of the plots shows the normalized differences with 
respect to the SHERPA prediction including spin correlations. 



and SHERPA with spin correlations on the one hand and 
of MCQNLO and SHERPA without spin correlations on the 
other hand show perfect agreement. Differences between 
the two spin correlation treatments may, thereby, easily 
reach up to 40%. 

The influence of spin correlations can also be seen in ob- 
servables based on two particle correlations. As two il- 
lustrative examples take the A<p and the AR distribution 
of the e + and the fi~ produced in the decay of the two 
W bosons. Again, the corresponding spectra, which have 
been exhibited in Figs. 1251 and 1261 differ significantly in 
shape depending on whether spin correlations are taken 
into account or not. 

The discussion of the impact of spin correlations is com- 
pleted by exploring the influence of the application of 
experimental cuts (cf. App. Q on the shape of certain 
spectra. It is clear that superimposing specific jet and 
lepton cuts strongly affects the event sample. Here, the 
cuts are mainly on the r\ and the pt of the leptons. In 
turn their distributions alter. The characteristics found 
for the cutfree case are not substantially changed by the 
applied cuts and by the renormalization of the spectra 
according to these given cuts indicated by the vertical 
lines in the Figs. 1231 and!24l More interestingly, however, 
these distributions drive alterations to secondary observ- 
ables. In the two-particle correlations mentioned before, 
the effects already present without applying cuts are en- 
forced. The slopes of the Atfi distributions increase, am- 
plifying the difference between both sets of predictions, 
the ones with and without spin correlations. The main 
change in the AR spectrum is an additional deviation 
between 0.2 (the cut) and 2.0, such that now the no- 
spin-correlation results are roughly 20% above the other 
ones. The case is different for the pseudo-rapidity dis- 
tribution of the W + boson. Without the application of 
cuts one starts off distributions that agree on the 10% 
level. This is severely changed by the introduction of 
the cuts, see the rightmost panel of Fig. [53 In contrast 
to the aforementioned two-particle correlations, here the 
predictions without spin correlations are well separated 
from the other ones only after the application of the cuts. 
As a last example, consider the transverse momentum 
distribution of the W + boson. Both types of predictions 



stemming from uncutted (left panel) and from samples 
analysed with cuts (middle panel) are pictured in Fig. 
1271 The inclusion of cuts apparently brings MCSNL0 and 
SHERPA including the full correlations into good agree- 
ment, but this clearly happened accidentally. 
To summarize, the examples shown here, clearly hint that 
the superposition of spin correlations (or their absence) 
together with cuts triggers sizeable effects in both types 
of observables, such that have already shown deviations 
in the absence of cuts and, more crucially, such that have 
not. In specific cases, such as the pt spectrum of the 
W + , this may possibly lead to misinterpretations of the 
results. 



V. CONCLUSION 

In this work, the merging procedure for multiparticle 
tree-level matrix elements and the parton shower imple- 
mented in SHERPA has been further validated; this time, 
the case of W pair production at the Fermilab Tevatron 
has been considered. First, it has been shown that the 
results obtained with SHERPA are widely independent of 
specific merging procedure details such as the choice of 
the merging scale and, for sufficiently inclusive observ- 
ables, the number of extra jets covered by the tree-level 
matrix elements. In addition, it has been shown that 
the specific form of the spectra produced by SHERPA is 
nearly independent - with deviations less than 20% - of 
the choice of the factorization scale and the renormaliza- 
tion scale. 

Having established the self-consistency of the SHERPA re- 
sults, they have been compared to those from an NLO 
calculation provided through MCFM. There, good agree- 
ment of the two codes has been found, again on the 20% 
level. Thus it is fair to state that the SHERPA results for 
the shapes are within theoretical errors consistent with 
an NLO calculation. The inclusion of the parton shower 
connected with specific scale choices in SHERPA, however, 
produces a surplus of QCD radiation with respect to the 
single parton emission in the real part of the NLO cor- 
rection in MCFM. 
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Finally, the results of SHERPA have been compared with 
those of other hadron-level event generators, namely with 
PYTHIA and MCONLO. In this comparison it turned out that 
SHERPA predicts a significant increase of QCD radiation 
with respect to the other two codes. For the pr spectra of 
jets accompanying the two W bosons, the differences are 
dramatic in the high-pT tails. In addition, the impact of 
spin correlations has been quantified. In the observables 
considered here, it reaches 20 . . . 50%. This may be even 
larger than the impact of higher order corrections. 
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APPENDIX A: INPUT PARAMETERS OF 
SHERPA 

All SHERPA studies have been carried out with the cteq61 
PDF set . The value of a s has been chosen according 
to the corresponding value of the selected PDF, namely 
a s = 0.118. The running of the strong coupling constant 
is determined by the corresponding two-loop equation, 
except for the SHERPA MCFM comparison. There an one- 
loop running has been employed for a s . Jets or initial 
partons are defined by gluons and all quarks but the top 
quark; this one is allowed to appear within the matrix 
elements only through the coupling of the W boson with 
the b quark. In the SHERPA MCFM comparison SHERPA 
runs, however are restricted to the light-flavour sector, 
i.e. the g, d, u, s, c sector. In the matrix element cal- 
culation the quarks are taken massless, only the shower 
will attach current masses to them. The shower cut-offs 
applied are 2 GeV and 1 GeV for the initial and the final 
state emissions, respectively. If explicitly stated a pri- 
mordial fcj_ Gaussian smearing has been employed with 
both, mean and standard deviation being equal to 0.8 
GeV. The Standard Model input parameters are: 

m w = 80.419 GeV , T w = 2.06 GeV, 

mz = 91.188 GeV , T z = 2.49 GeV, 

G M = 1.16639 x 10~ 5 GeV~ 2 , 

sin 2 6 W = 1 - m w /m%, 

a s = 0.118. (Al) 

The electromagnetic coupling is derived from the Fermi 
constant according to 

V2G, t M w sin 2 6w 
Oc em — . y-^^l 

The constant widths of the electroweak gauge bosons are 
introduced through the fixed-width scheme. The CKM 
matrix has been always taken diagonal. 



APPENDIX B: SETUPS FOR MCFM, MC@NLO 
AND PYTHIA 

MCFM 

The program version employed is MCFM v4.0. The 
process chosen is nproc=61. The investigations have 
been restricted to the d, u, s, c quark sector. The PDF 
set used is cteq61. The default scheme for defining the 
electroweak couplings has been used and their input 
values have been adjusted with the corresponding pa- 
rameter settings given for SHERPA. The renormalization 
scale and the factorization scale are fixed and set to 
Mfl = = Mw- 

MCQNLO 

The program version used is MC0NL0 2.31. The process 
number is taken as IPR0C=-12850, so that the under- 
lying event has not been taken into consideration. The 
two W boson decays into leptons are steered by the two 
M0DB0S variables being set to 2 and 3 for the first and 
the second choice, respectively. The lepton pairs have 
been generated in a mass window of 

M w - 40 T w < m 4u < M w + 40 T w . (Bl) 

Again, the cteq61 PDF set as provided by MCONLO's own 
PDF library is used. The weak gauge boson masses and 
widths are aligned to the settings used for the previous 
codes. All other parameters have been left unchanged 
with respect to their defaults. 

PYTHIA 

The PYTHIA version used is 6.214. The process 
pp -> W + W~ + X is selected through MSUB(25)=1. 
The specific decay modes of the two W's are picked by 
putting MDME(206,1)=2 and MDME(207, 1)=3, where all 
other available modes are set to zero. The possibility of 
parton shower emissions right up to the limit, which has 
been proven to be more convenient for jet production 
|55j . is achieved with MSTP(68)=2. This increases the 
IS shower start scale in PYTHIA to </I = 1960 GeV and 
accounts for a reasonably higher amount of hard QCD 
radiation. For all comparisons here, the underlying 
event is switched off, other parameters are left to their 
default. 



APPENDIX C: PHASE SPACE CUTS 

Two different analyses are used for the comparisons of 
the results obtained throughout this publication. A sim- 
ple analysis has been taken to verify the pure behaviour 
of the considered programs. For this case, only jets are 
analysed utilizing the Run II k±_ clustering algorithm de- 
fined in (5(| with a pseudo-cone size of R — 1. The jet 
transverse momentum has to be greater than 15 GeV. 
For more realistic experimental scenarios, an analysis ap- 
plying jet and lepton cuts has been availed. Then, the 
pseudo-cone size of the jet algorithm has been set to 
R = 0.7, and the jets have to fulfil the following con- 
straints on the pseudo-rapidity and the transverse mo- 
mentum, 

|?7> et |<2.0, #>15 GeV. (CI) 
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For the charged leptons the cuts on these observables are 
given by 

|r/ ep |<1.0, P £ p >20GeV, (C2) 

however, a cut on the missing transverse energy has not 
been introduced. There is a final selection criteria corre- 



sponding to the separation of the leptons from each other 
and from the jets, 

AR U > 0.2 , ARij > 0.4 . (C3) 
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